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Abstract. In this paper we obtain natural boundary conditions for a 
large class of variational problems with free boundary values. In compar- 
ison with the already existing examples, our framework displays complete 
freedom concerning the topology of Y — the manifold of dependent and 
independent variables underlying a given problem — as well as the order 
of its Lagrangian. Our result follows from the natural behavior, under 
boundary— friendly transformations, of an operator, similar to the Euler 
map, constructed in the context of relative horizontal forms on jet bundles 
(or Grassmann fibrations) over Y. Explicit examples of natural boundary 
conditions are obtained when Y is an (n-|- l)-dimensional domain in R"+^, 
and the Lagrangian is first-order (in particular, the hypersurface area). 
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Introduction 

Let y be a smooth (real) manifold of dimension n+1, with nonempty bound- 
ary dY. 

Definition 0.1. An n-dimensional submanifold L CY such that 

(1) L is connected, compact and oriented; 

(2) dL = LndY; 

(3) L is nowhere tangent to dY, 

is called admissible; the totality of such submanifolds is denoted by Ay- 

Introduce a local coordinate system {x, u) on Y, where x := {x^, . . . , x"). Let 
A = £d"a; be an r"^ order Lagrangian, i.e., let £ = £{x,u,ui) and / denote 
a multi-index of length < r. Suppose that, in such coordinates, an element 
L G Ay is the graph of a function u — u{x), defined on a connected and 
bounded domain C M": then the integral 

(0.1) Sx[L]:= J^^£(^x,u{x),^^d"x 

makes sense; it can also be given a coordinate-free formulation. 

Definition 0.2. The variational problem with free boundary values determined 
by the Lagrangian A on Y consists of finding the elements of Ay which are 
critical for S\. 

Indeed, if properly understood in a geometric framework, Sx is a real-valued 



function on Ay', the choice of the denomination is justified by (0.1): if L is 
allowed to vary within the class Ay, then the function u describing L is "free" 
to take any boundary value, as long as u maps dfl into dY. 

The main theoretical question addressed in this paper is the following: do the 
solutions to a variational problem with free boundary values fulfill some extra 
equation(s) besides the Euler- Lagrange equations^ A positive answer has already 
been given in [12l UHl |TT] , but without detailed proofs: Section |4] is devoted to 
review this result by adding the missing details. 

Sections [3}]4] deal with technical aspects of flag fibrations and relative C- 
spectral sequences, respectively: the reader not interested in theoretical consid- 



erations may skip them, and jump to Corollary 4.15 which summarizes their 
results. Section [l] explains the key used to obtain the main result (Section [5|, 
namely the natural behavior of the relative Euler map, under boundary-friendly 
transformations. As certainly know all who work in geometric variational cal- 
culus and cohomological theory of nonlinear PDEs, the Euler map is but a 
small feature of a general theory (comprising, e.g., conservation laws, Hclmoltz 
conditions, hamiltonian structures, recursion operators, etc.), which possesses a 
natural relative analog: we added Sections ISlB] just to give a ghmpse of it. 
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The applicative purpose of this paper is to present exphcit exaraples of natural 
boundary conditions. In the rather pedagogical Section[2] we review the classical 
analytic solution, given by van Brunt in a recent (2006) book to one of the 
simplest examples of variational problems with free boundary values. More 
involved examples are suggested by real-life circumstances, as, e.g., the problem 
of finding the equilibrium of a soap film freely sliding along the inner wall of an 



arbitrarily-shaped pipe, discussed in Subsection 5.2 



The geometric point of view is the backbone of this paper: besides allowing 
a transparent formulation of the main problem, it provides a key tool to obtain 
a solution. Analytic formulation (0.1) will be used whenever it is necessary 
to perform actual computations, as well as a source of valuable insights. For 
example, the Euler-Lagrange equations 

SC 

Su 

where ^ is the Euler-Lagrange derivative of C, are obtained by a well-known 



(0.2) 



0, 



manipulation of (0.1), under the assumption that the variations of u have a 

o 

compact support in fl: hence, a solution to the main problem should be a 
stronger condition than the Euler-Lagrange equations themselves. This clue was 
confirmed by the discovery of the relative Euler map |12|, reviewed in Section |4j 

1. Generalities on geometric calculus of variations 

The Euler map E appears, in one form or another, in all geometric frameworks 
for Variational Calculus that are based on the language of differential forms on 
jet spaces (called here Grassmann fibrations following the recent paper JT^ . of 
which we also adopt the notation). Building the Grassmann fibratiorFl G^^Y 
over Y is just a coordinate-free way to add new coordinates uj, with |/| < r, 
to the manifold Y, in such a way that A can be considered as an n-form on 
G^Y. In this new perspective, (0.1) can be rewritten without mentioning the 
local expression of L: indeed, since G^L = L, being dimi = n, the canonical 
inclusion 6/, : L C F is lifted to an immersion j'^L := G^z-l ■ L — > G^n^, which 



allows to pull any Lagrangian back to L. In other words, (0.1) reads 



(1.1) 



Sx-.LeAy 



fL*X e 



Passing from (0.1) to (1.1) is far from being a mere aesthetic exercise. It de- 



ploys powerful tools to attack the main problem: essentially, the possibility of 
using transformations which mix dependent and independent variables {x and 
M, respectively, in the above coordinate system). In Subsection 5.1 we show 



how a suitable change of coordinates can help avoid the lengthy computations 



^In Vinogradov and his school's approach, G^Y is denoted by J'^{Y,n), see I14| 
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proposed in Section [2j and how to obtain some useful formulae which, to the 
authors' opinion, would be very hard (though not impossible) to discover relying 
on pure analytic methods. 

The power of transformation methods descends from the natural character 
of the Euler map: in the principal geometric frameworks for Variational Calcu- 
lus (Krupka's variational sequences [1111], Anderson's variational bicomplex [T], 
and Vinogradov's C-spectral sequence [15]) the Euler map connects two spaces, 
say L{Y) and K{Y), containing, respectively, the Lagrangians and the Euler- 
Lagrange expressions for the variational problems on Y. We shall not go into 
the details, since a lot of excellent literature has been written on the subject; 
nonetheless, we stress that the natural character of the association Y i — > G^^Y, 
where r < oo, i.e., the canonical way to lift transformations of Y to the Grass- 
mann fibrations, makes the associations Y i — > L{Y), K{Y) natural as well. 
Indeed, L{Y) and K{Y) are usually defined as quotients of sub-complexes (or 
sub-sequences) of the de Rham complex of finite (or infinite-order) Grassmann 
fibrations, and as such they inherit the pull-back from differential forms. In 
other words, any diffeomorphism F : Y_ — > Y , determines a commutative dia- 
gram 

(1.2) LiY)-^KiY) 

F* f 

L{Y)^K{Y). 

If F is a wisely-chosen change of coordinates, then "the long way" from L{Y) to 
K{Y), i.e., F* o Ey o (F*)~^ may be more convenient concerning computations. 
But this is just a category-theoretic restatement of the well-known transfor- 
mation rule for the Euler-Lagrange equations, which was already known to E. 
Cartan: the purpose of this paper is to extend it to the class of variational prob- 
lems with free boundary values, where the Euler-Lagrange equations are sided 
by the so-called natural boundary conditions, or, equivalently, they are replaced 
by the relative Euler-Lagrange equations. 

Roughly speaking, the "relative" version of the Euler map arises because of 
the boundary dY. 

Definition 1.1. By abuse of notation]^ we shall put 

aC^r ip''-°)-\dY). 

Indeed, the canonical inclusion l : dG^^Y C G^^Y determines a differen- 
tial algebra epimorphism l* : fl{G'^^Y) — > r2(9GJ'jF) whose kernel ken* :— 
n{G';,Y, 9G;y) is, by definitionj^thc ideal of relative differential forms on G^-^Y . 



"^dG^Y is more like a prolongation, or lift, to O^Y, of the boundary dY. 
•^Such a construction is common in Differential Topology (see, e.g., [3]). 
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Much as L{Y), K{Y), and Ey are constructed out of (classes of) differential 
forms on G^Y and natural morphism connecting them, their "relative counter- 
parts", denoted by L{Y,dY), K{Y,dY), and Ey^, respectively, are built out of 
relative differential forms on G^Y. Details of this construction, carried out in 
the context of C-spectral sequences (meaning, in particular, r — oo), can be 
found in [H [TU] 

Section [4] explains why the relative Euler-Lagrange equations 

(1.3) E'^\X) = 

represent a solution to the main problem. More precisely, since K(Y, dY) iden- 
tifies with the direct sum K{Y)(BK{dY), the single equation ( |1.3[ ) captures two 
equations simultaneously, viz., the Euler-Lagrange equations 

(1.4) EyW^O, 

which involve n independent variables, and the natural boundary conditions 

(1.5) eU^)=0, 

where the number of independent variables involved is n — 1. Besides providing 
a common environment for such heterogeneous equations, the formalism of flag 
flbrations, introduced by the first author in II II, and reviewed in Section [s] 
allows to write down (1.5) in a workable wayrj 



The natural character of relative Euler map follows automatically from its 



very definition: in other words, the "relative" version of diagram (1.2), para- 
phrased by Lemma [l . 1 1 below, needs not to be proved. 

Lemma 1.1. Let F be boundary -friendly, i.e., F{dY) C dY. Then 

E'^^ ^ F* o Ef o {F*)-\ 

Lemma pTT] together with Corollary |4. 2. H will be employed in the last Section 
[5] to obtain new examples of natural boundary conditions. 

2. A MOTIVATING EXAMPLE 

Let n = 1, Y = [a,b] X R, and A = C{x,u,u')dx: in this case, functions can 
be identified with their graphs, and A := G°°{[a,b]) as a subset of Ay- Hence, 
up to a (noncritical) restriction of Ay to A, the boundary problem with free 
boundary values determined by A on Y, entails finding the functions u such thal[^ 

(2.1) lim ^#^^=0. 

||n-u||-!.0 \\u — u\\ 



"^This paper is based on the talk "A geometrical framework for Lagrangian theories which 
involve n and n — 1 independent variables simultaneously" delivered by the first author on 
August 24, 2012, within the conference "Variations on a Theme", dedicated to D. Krupka's 
seventieth birthday. 

^The norm can be either the L°° or the norm on C^{[a, b]). 
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In Chapter 7 of van Brunt's book [Ij, the above problem is modified by allowing 
u to be defined on a different interval than u. To fit this new setting, A must give 
up its linear structure and norm, namely A := '^xo<x-iC°° {[xq^Xi]); according, 
Y :— M? . Despite this, A keeps a rather obvious metric structure dji^. Moreover, 
with two real numbers Xq and Xi and a suitable function ^, one can construct a 
variation it of m, whose dyi-distance from u is controlled by a parameter e > 0. 
First, use Xq and Xi to define a new interval [io, ^i], where 

Xk^ Xk+ eXfe, fc = 0, 1, 

and suppose, without loss of generality, that xq = min{2;o,a;o} and 
xi = max{a;i, ail}. Then, use ^ € C°° {[xq,xi\) to construct the variation 

(2.2) u:=M* + e^, 

of u, where u* is the 2"^^ order polynomial extensiorj^of u to the interval [xq, ii]. 



u*{x) 

Now 



w, a; e [xq, Xi], 

u{xi) + {x ~ xi)u'{xi) + '^^'^^^ u"{xi), X € {xi,xi]. 



dA{u,u) := \\u~ u\\ + \[xq,uo) - (xo,uo)| + \{xi,ui) - (fi,ui) 



is a well-defined distance on A, which allows to adapt (2.1 1 to the case when the 



domain of definition of u can be altered: the norm \\u — u\\ has to be replaced 



by the distance d{u,u). Take the variation u (2.2 1, and compute 

(2.3) dAiu,u) < ||M-{t||+eXo + ||uo-wo||+eXi + ||Mi-ui|| < e(3C + Xo + Xi) 
Inequality ( |2.3[ ) shows that 

(2.4) dA{u,u)=o{e). 

In order to estimate the numerator in (2.1 1, comput^ 
(2.5) 

pxi pxi 

Sx[u] — Sx[u]= / C[u]dx ~ / £[M]da: 

Jxo Jxo 
/•xi+eXi i-xi 

— / £[M]da; — / £[u]dx 

J Xo+nXo Jxo 
pxi pXi+eXi pXo+eXo 

= / {£[u] - C[u])dx + / C[u]dx- / C[u]dx. 

J Xn J Xi J Xo 



Equality (2.5) shows that, with respect to the "fixed domain case" (2.1), the 



variation of Sx in u has two additional contributions due to the variations of the 



^To reduce the load of notations, we retain the same symbol u for the extension u* of u. 
"^It is convenient to write C[u] instead of C{x , u{x) , u' {x)) . 
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endpoints of the domain of u. The main advantage of the geometric approach 
presented in Subsection |5.1| later on, is that such a distinction between the 
variations of u and the variation of its domain, simply disappear. For the time 
being, (2.51 can be just rewritten in a more suggestive form 



where 



s s 

Sx[u] - Sx[u] = e^T^M 



{C[u] - C[u])dx 



SS\ SSx 



-I 









u] 


4 








J xa \ 



du dx du 



and 



£[M]da; = eXkC{xk, u{xk),u{xk)) + o(e^), A: = 0, 1, 



where we used the fact that 
d 
di 



Xk+t 



£[u]dx = C{x, u{xk), u{xk)) 

and C{x, u{xk), u'{xk)) — C{x, u{xk), u'{xk)) = o(e^), for fc = 0, 1. 

Now define real numbers Uo,Ui by eUk = u{xk) — u{xk), for fc = 0, 1, and 
compute 

eUk = {u + e^){xk) - u{xo) ^ u{xk) + eXku'{xk) + e^{xk) + o{e^) - u{xk) 
= e{Xku\xk)+^{xk)+o{e)). 

This shows that 

(2.6) ^{xk)^Uk-Xku'{xk)+o{e), fc = 0,l. 



In view of (2.6), (2.5) reads now 



dC 

S\[u] - Sx[u]=e < 



dC _ d dC 
du dx du' 



7 



^ (-1)^ -Xfe (^£-^.'1^^ {xk,u{xk),u\xk)) 



fc=0.1 



9£ a£ 

Uk^{xk,u{xk),u{xk)) +^{xk)Q;^{xk,u{xk),u{xk)) 



' dw 
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Equivalently, 
Sx[u] - Sx[u] 




{xk,u{xk),u'{xk)) \ +o(e^). 



Plugging the last expression into (2.1 1, and taking into account (2.4), wc finally 



see that u is a critical point for Sx if the above term in e vanishes for all variations 
of u, i.e., for all possible choices of ^, Xq and Xi. In particular, u must satisfy 
the (2"^^ order) Euler-Lagrange equations, 

SH 

— {x , u{x) , u' {x) , u" (x)) = 0, 
ou 

on its domain of definition, plus a (1^' order) natural boundary condition at the 
endpoints, 

dC 



(2.7) 



E 

fc=0,l 



{Uk~u'Xk)-^+XkC 



{xk,u{xk),u\xk)) = 0. 



Formula (2.7 1 is used in van Brunt's book to prove Theorem 2.1 below, which 



answers the main question for one of the simplest (though nontrivial) examples 
of a variational problem with free boundary values. 

Theorem 2.1 (Transversality conditions). Let F C be a closed and con- 
nected smooth domain, such that dY is the disjoint union of two curves 70 and 
71, and M.'^ \ Y is disconnected, and X be a 1*** order Lagrangian^ If an ele- 
ment u G Ay is a solution of the variational problem with free boundary values 
determined by \, then 

(1) u obeys the Euler-Lagrange equations on its domain of definition [xq, xi\; 

(2) u fulfills the following transversality conditions 
(2.8) 



dc7 



dx^ 



dcr 



c 



{xk,u[x^),u' {Xk)) = 0, 



0,1, 



□ 



where "fk[cr) = {x^^(a),y^^(a)) and 7fc(0) = {xk,u{xk)), fc = 0, 1. 
Proof. See [1], Chapter 7. 

In this Section we observed the lack of robustness of the functional-analytic 
approach: the slightest change of settings destroyed the norm on the class of 
admissible functions, and a (in many respects, unnatural) distance appeared in 
its place, which worked well only after some lengthy tricks. 



^Note that A^y is made precisely by all curves lying in Y , such that one endpoint belongs 
to 70 and the other one to 71, without being tangent to any of them. 
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3. Flag fibrations 



The main motivation for flag fibrations is that equations (1.4 1 and (1.51 in- 
volve n and n— 1 independent variables, respectively: merging them into a unique 
equation requires a new formahsm where the number of independent variables 
can take (at least) two values: n and n — 1. Recall the fundamental embedding 
G^nX ^ G\{GZr^Y). It allows to regard a point 9 G G^^Y as an n-dimensional 



tangent plane|jto G^^ ^Y, and an element of the fibered product 

(3.1) P:=G;r x^.-i^Gi_i(Grir) 

as a pair consisting of an n-dimensional and (n — l)-dimensional tangent plane 
to G^~^Y (at the same point). Define 

(3.2) K',n-l5^ := {(^n, On^i) G P | 0„ D 9„-l} . 

In many respects, the theory of flag fibrations parallels that of Grassmann fibra- 
tions; it is useful to review here some of its characteristic features. 

Theorem 3.1. Let r > 0. Then the following results hold: 

• F^^_iY is a smooth manifold, called the (r*^ order) flag fibration of 
Y (of signature {n,n — 1)): it is fibered over the base Y , as well as all 
lower-order flag fibrations, i.e., the manifolds F^^^^^Y, with s <r. 

• The flag fibration F^^_-^Y is naturally fibered over the corresponding 
(i.e., with the same order r and the same number of independent vari- 
ables n) Grassmann fibration G^^Y . 

• The image of F^ n-i^ under the canonical projection over Gl^_i{Gjj^^Y) 
is a smooth submanifold, naturally understood as 1^* order nonlinear par- 
tial differential equation on G^~^Y in n — 1 independent variables: the 
equation of involutive (n — l)-planes of G^~^Y. 

• The equation of involutive {n—l)-planes of G^Y projects naturally over 

n,n—l ■ 

• The infinite-order flag fibration F^^_-^Y , obtained as the inverse limit 
of finite-order flag fibrations, identifies with the equation of involutive 
(n — l)-planes ofG^Y, and, hence, it can be considered as an equation 
on G^r. 

• The infinite prolong atioi^^ (^7?°n-i^)(oo) of the 1^^ order differential 
equation F^^_iY , understood as a pro-finite leaf space)^ is naturally 
interpreted as a space of infinite-order Cauchy data. 

^Called integral element by Bryant&GrifEths 5 , or R-plane by Vinogradov and his school 

[HE]. 

^''See [7| for a definition of infinitely prolonged equations. 

^^In the sense of Vinogradov's "Secondary Calculus"; see, for instance, the introduction of 

nsi. 
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• The infinitely-prolonged equation (-P','f'„_ii^)(oo) fits into a double filtra- 
tion picture 



n ^ ^n-l-' • 

mimicking the similar diagram in the (linear) theory of flag manifolds. 

Proof See [U. □ 

Theorem 3.2 ([H], Theorem 9.1). Let L (resp., Y,) be a leaf (i.e., maximal 
integral submanifold with respect to the infinite-order contact distribution) of 
G^Y (resp., G^_iY). Then the following identifications 

p~\L) = G-i(L) 

(3.3) n'i(S) - J°°(iVs) 

hold, where Ns is a pro-finite vector bundle called the infinite-order normal 
bundle. 

Moreover, p and n are transverse one to another, in the sense that p^^{L) 
(resp., n~^(S)j maps non degenerately onto G'^_iY (resp., G'^Y). 



Equality ( |3.3[ ) is the less straightforward of the two, and plays a prominent 
role in the description of the relative Euler map, which will be introduced in the 
next section. 

4. Relative Euler operator and natural boundary conditions 

In order to clarify the relationship between relative cohomology and varia- 
tional problems with free boundary values, recall Definition |1.H and suppose 
that A = dXo, where Ao e fl"-^(G;;r) is such that 

Aolac^y = 0- 

Then reads 

(4.1) Sx[L]^ [ fL*X^ [ fL*dXo^ [ fL*Xo\oL = 0, 

Jl Jl JdL 



since dL C dY according to Definition 0.1 Indeed, j^L maps dL into dG^Y 



and, hence, the fact that Aq vanishes on the latter implies that its pull-back 
L* Xq vanishes on the former. 

Lemma 4.1. The action S\ on Ay is determined by the equivalence class of X 
modulo the subspace dn"-^{G';,Y,dG';,Y) ofn"{Gl,Y). 

Proof. A paraphrase of (14. 1|). □ 
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In order to simplify further analysis, we shall work, from now on, in the 
context of infinite Grassmann fibrations and C-spectral sequences; in particular, 
a Lagrangian will be a horizontal n-form on C^Y , 

where fth{G^Y) is the quotient differential algebra of fl{G^Y) with respect to 
the ideal of contact forms. An expert in bicomplexes or C-spectral sequences 
would say that the next corollary is the "horizontalization" of Lemma [4?!] above . 



Corollary 4.0.1. The action S\ on Ay is determined by the relative horizontal 
cohomology class 

niiG^Y) 



[XUeHj:{G^Y,dG:,Y) 
where d^ is the horizontal differential. 



dhni-'{G^Y,dG^Y)' 



Corollary |4.0.1| says precisely that H]l{G'^Y,dG''^Y) is the space L{Y,dY) 
mentioned in Section [l] The space K{Y,dY) can be obtained in a similar way, 
using relative forms, contact ideal, and cohomology: we shall rather use an 
approach based on total differential operators and Spencer cohomology, as in 
|12j . In the same cohomological framework it will also appear the relative Euler 



map Ey^, which allows to obtain the equation (1.3) out of the Lagrangian A. 



The aim of this section is to prove that (1.3) is indeed equivalent to the pair 



of equations (1.4)-(1.5) and, furthermore, that either the single equation (1.3) 



or the two coupled equations (1.4)-(1.5|, provide a (nontrivial) answer to the 
main question stated in the Introduction. The first result can be found in |12j . 



but its proof, which is a consequence of Theorem 3.2 was provided later in |11| . 
and it is a consequence of the following structural result, which dictates strong 
restrictions on the topology of the fibration dG^Y — > G^_i{dY). 

Corollary 4.0.2 ([E], Theorem 2). Consider G^_i{dY) as a submanifold of 
G^_i{Y) via the embedding dY C Y. Then 

dG^Y^Jf^{NG^_^(OY)) 



where N is the infinite- order normal bundle (see Theorem 3.2), and Jh means 
"horizontal jet bundle" 



It follows form Corollary 4.0.2 that the relative (i.e., constructed with relative 
forms) C-spectral sequence of dG^Y is particularly simple (i.e., one- linQ; 

m 

turn, this implies that K{Y,dY) splits into the sum K{Y) © K{dY) (the proof 



Roughly speaking, the analog of jet bundle where derivatives are replaced by total deriva- 
tives. See, e.g., 1141 1161 for more details. 
"'^'^See 1141 for the meaning of "one-line" 
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can be found in [TH])- Hence, equation (1.3) splits into two equations: (1.4) and 



(1.5) 



The second result has been stated in llj without proof, which is provided by 
Lemma 14.21 below. 



Remark 1. Proposition |4 . 1 1 below contains a general theoretical result concern- 
ing relative C-spectral sequences, so that there is no need to restrict ourselves to 
the case of one independent variable: in other words, we let Y to be of dimension 
n + m, where m is arbitrary, i.e., locally, to be fibered over an n~dimensional 
manifold X with m-dimensional fiber (when needed, such a fibration is called 
tt). Here we recall some terminology. 

VSym(y) is the module of vertical symmetries (denoted by >f in [HI?]) of 
the infinite-order contact distribution on G^Y, and V is the sub-algebra of 
differential operators generated by total derivatives (the C-differential opera- 
tors, according to [21 [7]). Suppose now we work in a local chart (in partic- 
ular, dY = {xn = 0} and tt is trivial): in this case, Dq^ denotes the pro- 
jection on the j"^ component of the free C°°(G^y)-module C°°(G^y)™ = 

VSym(y), and o'-f Di o D^/^ for all j = 1, . . . , m, and / multi-index of 
length n, i.e., / e N^. Moreover, V{C°°{G'^Y),ni{G'^Y)) identifies with 
V{G°°{G'^Y),C°°{G'^Y)) by means of the horizontal volume form d"a; , and 
V{G°°{dG^Y),ni-\dG^Y)) with V{C°°{dG^Y),C°°{dG^Y)) by means of 
the horizontal volume form d^~^x on dG^Y. Accordingly, the formally ad- 
joint modules (see [TB]) VSym^(F) and VSym^(9G'^y) are identified with the 
dual module of VSym(y) and VSym(9G^y), which are still free, with bases 

{Do'''}i=i,...,m and {£'o^^"'}j=i,...,„i,aeNo, respectively. 

Recall that Dj is the composition of total derivatives _D^\ °^^2 ° ' ' '"-D^" , with 
/ = (ii, . . . ,«„), and, by our own convention, the difference between the multi- 
index / and an integer a < i„ is the multi-index I ~ a :— (ii, . . . , in-i, in — o;). 

Proposition 4.1 (On the structure of K{Y,dY)). Let Y be as in Remark^ 
Then 

(4 2) K(Y dY) = V{VSymiYlKiG-Y)) 

^ ' ' ^ ' ' 6{V{VSym{Y),C^{G^Y))^ni-\G^Y,dG^Y))' 

where S is the Spencer differential. Moreover, the cohomology class of the cocycle 
□ = a^L'P^ e X'(VSym(r), flKG^Y)) is identified with the pair {E{D),E^{D)) 

e VSym^(r)® VSym^(aG5fy), where E{n) = Di{a^^)D^J\ and 

(4.3) 

m 
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Proof. By the definition of relative C-spectral sequences [12] , tfie space K{Y, dY) 
is tlie rt*^ cofiomology space of tfie subcomplex 



I?(VSym(y), C°°(G^r)) (g) nhiC^Y, dC^Y) 



of X>(VSym(y), n/i(G^y)). Expression (4.2) is a consequence of the fact that 
KiG^Y,dG^Y) equals ni{G^Y). In other words, K{Y,dY) has the same 
n-cocycles as K{Y), but fewer n-coboundaries, which explains why the n^^ 
cohomology of the subcomplex turns out to be quite larger than the cohomology 
of the entire complex: in turn, this explains the appearance of natural boundary 
conditions. For the sake of simplicity, we shall skip the index j. 

We now prove that the relative Spencer cohomology of □ is identified with 
(4.3 ). To this end, observe that the elements d"~^a;, for i = 1, . . . , n— 1, together 
with Xnd'^.-^x, form a basis for ni-\C°°{G^Y),dG^Y). Accordingly, ele- 
ments oiV{C°°{G^Y), ni-'^{C°°{G^Y), dG^Y)) can be obtained by summing 
up the elements D'^df-^a; and □"(g)a;„d;^-ia;, with e V{C°° {G^Y),C°° (G^Y)) 
for i = 1, . . . , n. 

We compute the Spencer differential ^(□^(g)d""^a;) = (-l)*-i(Di^on*)(g)d"a;, 
and observe that 



4(n"(/)) A ixr^d^-'x) + n'\f)d{xndii-'x) 

Di,^{a"{f))dXn A Xndn'^X + □"(/)dx„ A d'^^'^X 

(-l)"-i(x„Z?i„on" + n")(/)d"a;, 



for arbitrary / G G^{G^Y), whence 

(5(0" (g, x^dr'a;) = (-l)""'((x„i?i„ + 1) o □") (g, d"x 

= (-l)"-l(L>i„ OXnO □") d"x, 



since 1 — [Di^ , x„]. 

It follows that aoDi-oD is cohomologous to —Di^{a)oD for alH = 1, . . . , n—1, 
whereas a o Di^ o □ = (— (a) + o a) o □ is not generally cohomologous 
to ~Di^{a) o □, since Di^ o ao □ is not a coboundary, unless a factors through 

Xji- 
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Take now □ = a'Dj e V{C°° {G'^Y),ni), with / e N[J. Such an operator □ 
is cohomologous to the operator 

□' = (-i)'""*"^/-»„(«')°^;: 

= {-l)\'\-'-{~D^,SDI-^Aa')) + i^.„ o (a^)) o I?;--! 

: in iterations 

= {-l)\'\Di{a') + (-1)1^1-'" i?,„ o i?/_,„(a^) o D'--^ 
i.e., to 

W = {-lf\Dj{a')+ {-l)\'\-^D,^oDj^^{a')oD2:\ 

0<Q<i„ 

which turns into a coboundary if and only if the function Dj{a^) is zero and all 
the functions Dj^aia^) factor through x„, i.e., they vanish on dG^Y. 

This means that the cohomology class [□] is uniquely determined by 
(-l)l^l£i/(a'f) e C°°{G^Y), i.e., by E{a) and by the set of functions 
l-l)\i\-°'Di_a{a^)\OG^Y e C°°{dG^Y), with a = l,...,i„. Notice that the 
latter ones can be rewritten as Dj-i^{D^^{a')^QQa^Y), with a = 0, . . . , i„ — 1, 
because the multi-index I — in belongs actually to Nq~^ and the first n—1 total 
derivatives are tangent to dG^Y: hence the last set of functions represent the 
coordinates of the element E^{\I\), according to (4.3). □ 

Proposition 4.2 (On natural boundary conditions). Let L be locally given by 
the graph of a function u ~ u{x), defined on a compact and connected subset 
n C M", such that dfl has equation x" = 0. If L is critical for S\, then the 
following equations hold on dfl: 



dC 



r\\I\-i, 



an 



= 0, a = 0, . . . , r — 1. 



Proof. Let r be the order of A and I — (ii, . . . , i„). Put 

dui du^ii^i2 
From the well-known formula of elementary calculus 

(4.4) /gW = i^J^'^a^'^'^'^^ + (-l)7^^^i 
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it follows that 



' iiH \-in<r 



(4.5) / Yl ^'''-'"V^^^^ 
(4.6) 

(4.7) +((-1)*^ 



E 

ilH hiji'^r 



d 



\si=0 



d"cc 



Then, applying again (4.4 1 to the term (4.7), we obtain 



(4.5) =(^4.6 



(4.8) 
(4.9) 



(-1) 



E 

i'l H \-in 



d 

da;2 



'12-1 



*1 . 



,*1 



{-ly 



E(-ir 

Jii£iii2---i„ 



d«2 
d^ 



da;!^ 



d"a; 



Again, by (4.4), we develop term (4.9): 



(4.10) (4.5) ^ (|4.6|) + (|4.8 



E 



iiH h-'in^'" 



(4.11) 



(4.12) 



1 da;„ 
(-ir 



' ^E(-ir"S:((-ir 



\s„=0 



dx^" 



+(-ir 



(-1) 



dx\' 



77 \d''x. 



Since 957 = {a;„ = 0}, all terms appearing on line (4.101 disappear, being of the 
form 



/ ^(F)d"a;= / d{Fd'^-^x)= [ Fd'^-^x] 
J ft (^^i Jn Jon 



0, Vi = l,2,...,n- 1. 



On the other hand, (4.12) is the Euler-Lagrange; it remains just (|4.1l|), i.e 

dl^l-' 



(4.13) 



ion 



E E(-i)'^'- 

|/|<r s^-O 



■in-\-Sn 



dx^-*" 



dXr. 



dr^a;. 



an 
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Since (4.13) has to vanish for all variations rj, all equations (4.16) must be 



satisfied. 



□ 



Lemma 4.2. Let L G Ay be a solution to the variational problem with free 
boundary values determined by X on Y . Then equation (1.4) holds on L and 



equation (1.5) holds on dL. 



Before providing a proof, it is convenient to cast a bridge between the ap- 
proach based on total differential operators to the space K{Y), sketched in 
Reniark[l] and a perhaps more familiar one, based on "1-contact, n-horizontal" 
(n + l)~forms, or forms "of type (l,n)". Namely, (1.4) can be written down as 



(4.14) 



Ey{X) = — cjAd"a;, 
ou 



where w is the zero-order contact form, and uj A d"a; plays the role of the gen- 
erator Dq (see in Remark[T]) of the module K{Y). Equation (4.14) clarifies the 
above sentence "(1.4) holds on L": it means that (4.14), pulled back to L via 
j°°L, vanishes. 

Similarly, the results contained in Corollary |4.0.2| and Proposition |4. 1| give a 
solid basis to the sentence "holds on 9L" , since 



(4.15) 

where now the 



1 



w"'s are the zero-order contact forms on Jj^{NQ=o ^(gy)), and 

w" A d^-iX plays the role of the generators Dq" , where j ~ 1 (see in Remark 
[T]), of the module K{dG^Y). 

According, EyaiX) can be obtained as the coefficient of -DJ'", in (4.3), where 



j = 1, and — Again, the sentence "(1.5) holds on dL" means that (4.15) 



pulled back to dL via joodL, vanishes (see also Theorem 11.1 in |11)). 

Proof. The first fact is obvious: if L a solution of a variational problem with 
free boundary values, then L \ dL is a solution of the Euler-Lagrange equation 
determined by the same Lagrangian A, i.e., equation (1.4) holds. 



We stress that, in order to prove the second fact, it is necessary to have the 

Namely, 



result on the structure of equation (1.51 provided by Corollary 4.0.2 



(1.5) is localizable, in the sense that its left-hand side belong to a module of 



sections, and hence it vanishes locally if and only if it vanishes globally. Then, 
we can choose a coordinate system {x, u) such that L is the graph of a function 
u = u{x) on Q, and dY — {x^ = 0}. Since L is critical, all equations (4.16) 
must hold true on 917; on the other hand, the above discussions showed that 



equations (4.161 are nothing but E 



d 



0: hence, (4.15) vanishes, i.e., (1.5) 



must be valid on dL. 



□ 
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For readers not interested in theoretical details, we collect the main result of 
the last two sections into a convenient (though redundant) Corollary. 

Corollary 4.2.1 (A solution of the main problem). Let L £ Ay be a solution 
to the variational problem with free boundary values determined by A on Y. 
Then, the natural boundary conditions EyiX) = are satisfied on dL. In local 
coordinates, Eyi^) ~ Ey^aW^" ^ ^n-i^r where: 

i) the uj" 's are the zero-order contact forms on the infinite jet of a suit- 
able pro-finite vector bundle over dY which arises in the theory of flag 
fibrations over Y ; 

ii) if a is less than the order of the Lagrangian \, and L is locally the graph 
of a function u = u{x) on il <Z M", the component Ey^[X) is given by 



5. Applications 



Together, Lemma |l.l| and Corollary 4.2.1 provide a powerful tool for writing 



down concrete examples of natural boundary conditions. Computations pre- 
sented in this section will be simplified by some "tricks" based on multi-linear 
algebra and total differentials (Remarks [2] and [s] below). 

Remark 2 (Top differential forms). A brute-force attempt to change variables 
in a multi-dimensional integral may lead to a meaningless formula 

(5.1) j f{x)d-x = I f{xit))'^dH. 

Nonetheless, since the C°°(M")-module f2"(IR") is freely generated by d"t, any 
n-form can be identified with a function. In particular, d"t identifies with 1, 
and d"a; with the Jacobian of the change of variables x — x{t), thus recovering 



the meaning of (5.1 ). From now on, all n-forms will be identified with functions: 
hence, an expression like S(cli), where 2 is a vector field and w an n-form, is not 
the Lie derivative of uj, but the function S(/), where / is uniquely defined by 
w = /d"t. 



Remark 3 (Total differentials). Formula (5.1) can be adapted to variational 



integrals, just by replacing differentials by total differentials, namely 

f{x,u{x),Ui{x), . . .,Unix))d"x = 



d X 

f{x(t,y),u(t,y),ui{t,y,yi, . . .,yn),-- . , u„(t, y, yi, . . . , ?/„))^j^-d"t 

d t 
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where now d"x — dx^ A • • • A dx" (i.e., the operator dh used in Section|4|. Recall 
that 

dx = Dii{x)df, X — x(t) 

where Dfi ~ dfi + Uidu is the total derivative operator with respect to t^. 
In this context, horizontal n-forms on Gl^Y, i.e., the space with coordinates 
(t, y,yi . . . , Un), are identified with functions on the same space, via the horizon- 
tal volume form d t. Accordingly, d a; is the "total Jacobian" associated with 
the change of variables {x,u) = {x{t,y),u{t,y)). 

5.1. A 1*** order, one— dimensional example. Consider again the variational 
problem with free boundary values of Theorem 2.1 Section[2j Let C = C{t, y, y') 
be its Lagrangian, and recall that dY is the disjoint union of two curves in the 
(t, y)-plane. Then, if 7(ct) = {t^{a),y^{a)) is one of them, a critical point 
y = y{x) for Sx must fulfill the natural boundary condition 



(5.2) 

where (t, y{t)) 



y'-rio) 



dc 



dy' 



C 



{t,y{t),y'{t))=0 



7(0) (see (2.8 1). 



Equation (5.2) can be obtained in a transparent geometrical way, without 
introducing ad hoc metrics on the set A. Just use a change of coordinates 

X = x{t,y) 
u = uit.y) 

which "rectifies" the curve 7, i.e., such that i^*(7) is, for instance, the w-axis of 
the (cc, u)-plane. Then, lift to a contact transformation (t,y,y') 1 — > {x,u,u') 
of G\Y, where 



(5.3) 



ut + yuy 
xt + y'xy ' 



It is a simple exercise to get (5.3 ) (see, for instance, 2J, Section 1.2); nevertheless, 
in view of the next generalization, we prefer to justify it, by using the total 
differential operator d. Recall that 



(5.4) 



df^Dt{f)dt, 



for / = f{t,y,y'), with Dt — dt + y' dy being the total derivative operator in 
t. As announced in Remark [3j we shall identify horizontal one-forms on the 
(t, y, ?/')-space with functions: hence, (5.3) above can be written as 



(5.5) 



du 
dx 



Dt{u)dt 
Dt{x)dt 
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It is worth noticing that the inverse transformation ^ is the same as (5.5) 
(5.6) 



y 



dy Dx{y)dx _ Vx + u'yu 



dt DAt)dx tx+u'U' 



where now the total derivative operator is taken with respect to x. It foUows 
that 



(5.7) 



du' 



d_ Z D Ay) 

du' \Dx{t) 



Finally, since dt = dt (see (5.4)), the Lagrangian A reads 

A = £{t, y, y')dt = C{t, y, y')Dx{t)Ax 

in the (x, m, M')-space, where Dx{t) plays the role of "total Jacobian" (Remark 
In other words, (F*)~^(A) = Cdx, where C = C{x,u,u') is given by 



(5.8) C ^ {F*)-\C)Dx{t) = C{t{x,u),y{x,y),y' {x,u,u'))Dx{t){x,u,u'). 

Now everything is ready. {F*)^^{X) determines a variational problem with free 
boundary values on F{Y) and dF{Y) — F{dY), by construction, consists of 
two curves, one of which is the w-axis: by Corollary |4.2.H on such an axis the 
natural boundary conditions take a particularly simple form 



dC 

du 



^(0,u(0),w'(0))=0. 



It remains to express (5.9) in terms of the coordinates {t,y,y'), i.e., to apply 
Lemma [LT] We compute 



du' 



du 



- dy'du' U.wy ^ 



du' 



dCDx{t)^{Dx{y))-Dx{y)^{Dx{t)) , BD^it) 



dy' 
dC 



Dx{t) 



+ C 

c 



du' 

dDxlt) 



(5.9) 



5^ dC 

dy 



du' 
dDx{t) 
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It suffices to observe tliat 



(5.10) 
(5.11) 



du' 
du' 



d{U + u'tu) 

du' 
djVx + u'yu) 

du' 



Substituting (5.10) and (5.10) into (5.9), one gets 



dC 

dy 



7 {Uu - tuy') + Ctu = 0, 



which, evaluated at (0, u(0), u'(0)), returns (5.2), since, by the choice of f , 

t;(o) = i„(7(o)), 

y;(0) = y„(7(0)). 

5.2. A 1*** order, multi dimensional example. Now we pass to an n-dimen- 
-sional example: as we shall see, the geometric methods used before generalize 
effortlessly to this case; an analogous generalization of the methods used in 
Section [2j i.e., defining a metric structure on the space of all functions defined 
on a compact and connected subset of M", would introduce a lot of technical 
difficulties, obscuring the simple solution of the problem. 

Theorem 5.1 (Generalized transversality conditions). Let Y be a closed smooth 
domain, with nonempty (smooth) boundary, of M"^^ — {t,y), with 
t ~ {t^, . . . ,t"), and X be a 1^*' order Lagrangian, locally given by 
C = C{t,y,yi, . . . ,yn). If y — y{t) is a critical point for S\, then, for any 
point 6 g dY , the normal vector vg to the hypersurface dY must be orthogonal 
to H{t, y,y^(t), . . . , y"-(t)), where H is the M"~^^ -valued function on G\Y given 

by 

„ fdc dc ^ dc' 

V oyi dyn 



and e = {t,y{t)). 

F 

Proof. Choose a change of coordinates (t, y) \ — 

X = x{t,y), 
u = u{t,y), 



dy. 
ix,u), 



such that F{dY) has equation a;" — 0. In analogy with (5.6), 
(5.12) y, 
where 



dy 
dp 



~i — :n— 1 , 

dy Ad- t 

— . —71—1 

dP A d.. t 



dt^ A • • • A df-' AdyA df+' A • • • di" 
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and we use again the convention that horizontal n-forms are identified with func- 
tions via the (horizontal) volume form dx" introduced in Remark [sj Developing 
all total differentials appearing in (5.12), one recovers the familiar formula for 



the lifting of F, as it can be found, e.g. in 2], Section 1.2. 



In analogy with (5.7) 



dyi _ _d_ f 



i = l,2,. 



. ,n. 



Finally, analogously to (5.8), we obtain (F*) ^(A) = Cd"x, where 
C = C{x, u,ui, . . . , Un) is given by 



(5.13) 



C^{F*)-\C)d t 



where, as explained by Remark |3j c? t = =k^ is just an unconventional way to 
write down the "total Jacobian" of F. 

Again, {F*)~^(X) determines a variational problem with free boundary values 
on F{Y) and dF(Y) — F{dY), by construction, is the hyperplane x" = 0: by 
Corollary |4. 2. 1[ on such a hyperplane the natural boundary conditions read 



,0,m(x\ 



,0), 



.,x 



,0))-0. 



Finally, let us write down (5.14) in terms of the coordinates (i,y,?/i. 
We compute 

OUn V / 



,yn)- 



du 



dUn dUn 

% dUn Xd^tJ dUn 

^ dC d"t^iu;,)-u;,^{d''t) ^ ^ d{Tt) 

dUn 



dC f d 



d t 



%i V^"" d tdun 



C 



d{d t) 

dUn 



(5.15) 
We denote 



5^ d£ j duji 
dyi \ dun 



y 



d{d"t) \ ^ 

dUn j 

dujn d{d'^t) 
' du„ ' du„ 



L 



d{d"t) 

dUn 



Then, (5.15) coincides with v ■ H. 
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It remains to show that v is indeed the normal vector to F{dY). To this end, 
it is convenient to pass to the determinantal notation for n-forms. Namely, 



Vxi +uiyu 

+ U2Vu 
yx„ + UnVu 



tx,+ 



contains w„ only in the last line: hence, 



dijjj 
du„ 



Vxi 

yx2 



U2yn 



4-n I 4-n 



tx, + Ultl 
tL + 



yx,^- 



t. 



Subtracting from the j"^ row the n"^ row multiplied by u.j , for all j = 1, , 
the determinant does not change, i.e.. 



(5.16) 



Similarly, 



(5.17) 



dUn 



tL 
tL 



yx2 
yu 



d{d t) 

dUn 



xi 



tL-. 
tl 



Xn-l 



j-n 

^x„-i 
''11 



Observe that (5.16) and (5.17) are the multi-dimensional analogues of (5.10) 



and (5.11), respectively. In other words, v is composed of the n x n minors 



(with sign) of the n x (n + 1) matrix 

{Ti,T2, . . . , T„_i, T)*, 

where the n vectors 

Tl = (t^i,y^.), i = l,...,n-l, 

T i^u: yu\ 

form a basis for the tangent space of F{dY) — {xn — 0}, i.e., v = Ti x ■ ■ ■ x 

Tn-l X T. □ 
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Notice that (5.15) is formally the same as (5.9): the synthetic language of 



multi-linear algebra allowed to handle all the "total Jacobian" determinants 
involved in the proof, without any additional difficulty as compared with the 
one-dimensional case. 

5.3. The soap film. We generalize now a classical example that can be found in 
Giaquinta and Hildebrandt's book ^ (Section 2.4). Namely, in the hypotheses 



of Theorem 5.1 above, suppose that A is the (hypersurface) area Lagrangian, 
i.e., locally. 



Then 



H = 



1 



1 



(yi,2/2,---,2/ri,-l) 



is precisely the unit normal vector to the surface y = y{t). This proves the last 
result of this paper. 

Corollary 5.1.1. Let Y C M"+^ be a closed smooth domain with smooth 
nonempty boundary. If a hypersurface L G Ay is a solution of the variational 
problem with free boundary values determined by the area functional, then L 
must intersect orthogonally dY everywhere. 



In particular. Corollary |5 . 1 . 1 1 shows that a soap film, whose boundary is con- 
strained to slide over the inner surface of a fixed domain (e.g., a pipe of arbitrary 
shape), tends toward a position of equilibrium where it forms a right angle with 
the walls of the container (besides, of course, possessing zero mean curvature). 

Remark 4. If n = 2 and Y = DxM.'^, where D C M is diffeoniorphic to a closed 
disk, then a surface L from Corollary |5.1.1| above is forced to be the graph of 
a constant function y = y(ti,t2). Indeed, since 1" is a surface with zero mean 
curvature, its maximum is attained on dD. Hence, there exists a point 6 e dL, 
such that dL = L O dY has negative curvature. But L hits dY orthogonally in 
9, thus, along the normal direction to dL, the surface L must possess positive 
curvature, i.e., there must exist a point 9' € L, in a neighborhood of 9, such that 
the t/-component of 9' is greater than the y-component of 9, thus contradicting 
the fact that 9 corresponds to a maximum. It follows that L must be the graph 
of a constant function. It would be nice to generalize this simple observation to 
multi-dimensional cases. 
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